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Abstract 

We present a novel analytic calculation of the Haar power spectra, and various higher order 
spectra, of mean field avalanche models. We also compute these spectra from a simulation of the 



' results and obtain novel exponents. 
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Introduction: Certain disordered systems crackle, that is when such a system is slowly 
driven it responds by discrete jumps or avalanches of a broad range of sizes This 
crackling is often the result of non-equilibrium collective transport in random media, where 
the interactions are strong enough that thermal effects are negligible and the system has 
many interacting degrees of freedom. Such systems include: charge density waves, vortices in 
type II superconductors, crack propagation, earthquakes, and Barkhausen noise in magnets. 
In mean field theory the avalanche dynamics of the systems mentioned above are described 
ry Possion statistics (given by Eq. ((T])) at the critical point or critical depinning transition 

m. 

In our novel mean field analysis we determine the following spectral functions: Haar 
power spectra, 1.5 spectra, and second spectra for systems that have avalanche dynamics 
given by Eq. (^. Haar power spectra allow us to obtain a power series in time, H{t,fi), 
needed to calculate higher order spectra. Higher order spectra give valuable information 
about the avalanche dynamics in Barkhausen noise not accessible through ordinary power 
spectra QS. Also, higher order spectra have bee,r used to obtahr crucial informatiorr about 
a varietyof diverse systems such as: metastable states in vortex flow p], natural auditory 
signals ^ conductance- noise power fluctuations [8J, non-Gaussian 1// noise 0], and spin 

[lasses IC L While much experimental work has been done studying higher order spectra 
ifly,[aB,y,[io 



range of systems 



we present a rigorous mean field treatment that is applicable to a broad 
j|. This analysis is one of a small number of theoretical calculations in 
higher order noise statistics. 

We also compare our general results from mean field theory to a mean field simulation 
of the T = random field Ising model (RFIM) and an analysis of Barkhausen noise 
obtained from a simulation of the T = infinite range RFIM in ci = 3. We also find novel 
exponents from our analysis. 

Mean Field RFIM: The T = mean field RFIM consists of an array of spins (sj = ±1), 
which may point up (sj = +1) or down (sj = —1). Spins are coupled to all other spins 
(through a ferromagnetic exchange interaction J), and to an external field II{t) which is 
increased adiabatically slowly. To model dirt in the material, we assign a random field, hi, to 
each spin, chosen from a distribution P{hi) = exp{—hj/2R'^)/\/27rR, where R, the disorder, 
determines the width of the Gaussian probability distribution and therefore gives a measure 
of the amount of quenched disorder for the system. The Hamiltonian for the system at a time 
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t is given by: H = — J2{JM + H(t) + hi)si, where M = is the magnetization of the 

i j 

system. Initially, H{—oo) = — oo and all the spins are pointing down. Each spin is always 
aligned with its local effective field h1^^ = JM + H{t) + hi. 

RFIM with Infinite Range Forces in 3D: The T = Infinite Range RFIM (IRM) consists 
of a 3D lattice of spins, with a Hamiltonian at a time t is given by: H = —JsiSj — 

<ij> 

J2i{H(t) + hi — JinfM)si, where Jj„/ > is the strength of the infinite range demagnetizing 
field, and (ij) stands for nearest neighbor pairs of spins. The local effective field is given by 
hi" = JE<ij> Sj + Hit) + hi- JinfM. 

The addition of a weak Jj„j ~ ^ to the traditional RFIM causes the system to exhibit 
self-organized criticality (SOC) This means that as H is increased the model 

always operates at the critical depinning point, and no parameters need to be tuned to 
exhibit critical scaling behavior (except ^ 0). We limit our analysis to a window of 
H values where the slope of M{H) is constant and the system displays front propagation 
behavior. Details of the simulation are given elsewhere (l5| . 

Avalanche Dynamics in the T = RFIM: The external field II{t) is adiabatically slowly 



increased 



to flip 



rom — oo until the local field, hi" , of any spin Si changes sign, causing the spin 



16|. It takes some microscopic time At for a spin to flip. The spin flip changes 



the local field of the coupled spins and may cause them to flip as well, etc. This avalanche 
process continues until no more spin flips are triggered. Each step of the avalanche, that 
is each At, in which a set of spins simultaneously flip, is called a shell. The number of 
spins that flip in a shell is directly proportional to the voltage V{t) during the interval At 
that an experimentalist would measure in a pick-up coil wound around the sample. In our 
simulations we denote the number of spins flipped in a shell at a time t by nj(= V{t)). 
The first shell of an avalanche (one spin flip) is triggered by the external field II{t), while 
each subsequent shell within the avalanche is triggered only by the previous shell, since 
H{t) is kept constant while the avalanche is propagating. H{t) is only increased when the 
current avalanche has Stopped, and is increased only until the next avalanche is triggered 
(i.e. ^ 0). The number of shells in an avalanche times At defines the pulse duration, 
T, or the time it took for the entire avalanche to flip. The time series of Ut values for many 
successive avalanches creates a Barkhausen train. 

Poisson Distribution: In order to calculate the Haar power spectrum and higher order 
spectra we first need the probability distribution for an (infinite) avalanche at criticality in 
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the class of mean field avalanche models we are studying j3[ , which is given by the following 
Poisson distribution 



An avalanche begins with a discrete jump of unit magnitude, so we have = 1, in 
the context of Barkhausen noise each avalanche begins with one spin flip. Now (...) = 
J2 ...P{ni, ...rioo) represents the average over Eq. (^. Using Eq. we determine the 

{ni,...,noo} 

2-pt, 3-pt, and 4-pt time-time correlation functions, that allow us to analytically calculate 
the Haar spectra and higher order spectra. Details are left for a long paper. 

Haar Series and Haar Transform: The Haar transform, H{t, /i), a simple wavelet trans- 
formation, is the square wave equivalent of the Fourier transform; where H{t, fi) is the 
absolute square of the time integral over a period (of duration j^) of the square wave times 
a section of the train centered around t |4]. The Haar transform affords us the needed time 
resolution at the expense of frequency resolution that is extraneous to our purpose, since 
we are interested in studying how the power contribution at a given frequency /i changes 
in time. 

We first determine the Haar power series for a single large avalanche of duration T (in 
the limit that T oo) using the 2-pt correlation function. We obtain the following exact 
result: H(t, fi) = -^(2^^^ + 4-) Moreover, to find the Haar power we sum our result over 
the whole avalanche to obtain: 



This result differs from the Fourier power, Spifi) = ji [2], by an additive constant and 
a constant factor. This is an artifact of doing a discrete transform. We sum from t = 1 to 
T/i (not to T) since we are summing over Haar wavelets of duration j^. We reaffirm the 
above result (Eq. through a mean field simulation. 

1.5 Spectra and Second Spectra: The 1.5 spectra, and second spectra are defined below: 



oo 

P(no = 1, ni, ris, ^oo) = r TT 

eni! fj^ 
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where v(t, /i) is the sum of Ui around time t over a duration j^, and /2 is the frequency 
conjugate to t. Also, Ft{...} is the discrete Fourier transform (FT) with respect to t, and 

To calculate Eq. (jHlEl) we first write the product of FTs as the FT of a convolution. 
This leaves us with the following quantities: {v{t, fi)H(t + 9, /i)), and {H(t, fi)H{t + 9, fi)) 
where 9 is the convolution variable. These quantities may then be rewritten as sums of 3-pt 
or 4-pt correlation functions, and subsequently evaluated. We obtain the following results: 



ri.5(/i) = E(^^(t,/iW,/i))^^ (5) 
t=i Jl 

T,Ui)^Y.mJi)H{tJi))^^, (6) 
t=l J I 

Re{SiM2ji)} = ^ + r^i2ifi) (7) 

^2(/2,/i) = || + rf)(A) (8) 

h 

where V1.5, V2, Q1.5, Q2 are given in Table IJ for mean field theory, the mean field simu- 
lation, and the IRM. Also, ^1,5, A2, -B1.5, and B2 are non-universal constants. Plots of Eqs. 
(I3IH1) are given in Figs. ©-(EI)- 

Now ri/2(/i) and r2(/i) are independent of /2, and result from the ^ = evaluation of 
(f (t, fi)H{t + 9, /i)), and {H{t, fi)H{t + 9, /i)^ (see Fig. [l}. These /a independent terms 
are referred to as the DC background terms [9]. The functions r[^''(/i) and r2^^(/i) are 
ri.5(/i) and r2(/i) normalized by J2t H(t, f\) and (X)* H{t, /i))^, respectively. Consideration 
of Irn{Si,5{f2, fi)} is left for a long paper. 

We have left Eq. (jSHHI) in terms of general scaling forms since we ascertain that our mean 
field and IRM simulations obey the same scaling forms only with different exponents and 
non-universal constants. 

Mean Field Simulation: We perform 300 runs of a simulation of the mean field RFIM. 
We collect data taken from H G [0, 0.00125] at i? = 0.79788 {Rc = 0.79788456) in systems 
with = 15 X 10® spins, and J = 1. The results of our simulation agree with the scaling 
forms given in Eq. (jSHH)) with exponents given by Table IH See Fig. 

Infinite Range Model Simulation: We perform 60 runs of a 3D simulation of the infinite 
range RFIM. The data was taken from H G [1.25, 1.88] (from the slanted part of the hys- 
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^1.5 


V2 


Ql.5 


Q2 


MFT 


2 


2 


3 


5 


MF Sim. 


1.93 ±0.10 


1.92 ±0.12 


2.95 ±0.10 


4.93 ±0.12 


IRM {d = 3) 


1.80 ±0.07 


1.80 ±0.05 


2.73 ± 0.06 


4.46 ± 0.07 



TABLE I: We present the values of the exponents: Vi,^, V2, Q1.5, and Q2 given in Eqs. for 
mean field theory (MFT), mean field simulation (MF Sim.), and the infinite range model (IRM) 
for d = 3. While the exponents values for MFT were determined analytically, the exponents for 
the MF Sim. and IRM were determined through a non-linear curve fitting of simulation data. 

teresis loop) at R = 2.2 in system with = 400^ spins, J^f = j^, and J = 1. Again, the 
results agree with Eq. (jSHHl) with exponents given in Table ID Refer to Fig. [T]|ni 

Finite Size Effects: We check finite size effects in our mean field and IRM simulations by 
looking at higher order spectra for various system sizes. We find that for smaller system sizes 
the high frequency scaling (and flattening due to the DC background term) is unchanged 
for second spectra and 1.5 spectra. However, at low frequency the scaling regime of the 
second spectra and 1.5 spectra flattens (in MF Sim.) or rolls over (in IRM) at frequency 
fi ^ 7r^, where T^ax is the maximum avalanche duration. Tmax is system size dependent 

J max 

since Tmax ~ L^, where L = N^^'^ (d = dimension of system). For a N = 400'^ system we 
flnd that T^ax ^ 4300At. 

Discussion: We flrst notice that the exponents for our MF simulation are systematically 
smaller (by an amount of 1% to 4%) than the MFT exponents, since the MF simulation 
results were obtained from a train of flnite avalanches while our MFT calculation was for a 
single inflnite avalanche. The intermittency introduced by a train of flnite avalanches appar- 
ently lowers the magnitude of the mean fleld exponents since the intermittency effectively 
adds white inter-avalanche noise to the intra-avalanche noise seen in the MF calculation. 

For the DC components {6 = 0) of the higher order spectra (see Fig. ^ we flnd excellent 
agreement between MFT apd the MF simulation. In fact the exponents and Q,, are 
not novel; with the help of llj we ascertain the following exponent relations: Q2 = — 2 
and Qi,5 = ± 1. Using the values given in 0| for r and Xjavz we flnd exact agreement 
in MFT: Q2 = ^ and Q1.5 = 3, and in the IRM in = 3 we flnd Q2 = 4.40 ± 0.10 and 
Qi s = 2.72 ± 0.03, in close agreement with the table above. 
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In Fig. (j21ISI) we present the 1.5 spectra and the second spectra. We determine the high 
frequency power law exponent for the 1.5 spectra and second spectra given by the exponents 
Vi.5 and V2 respectively (see Table H)). Upon subtracting the DC background we find that 
the second spectra, for different fi values, collapses (see Inset of Fig. Since the DC 
background term is small for the 1.5 spectra, we find a collapse without subtracting the DC 
background. 

The weak dependence of 1.5 spectra on /i and the strong fall off of the second spectra 
at high frequency in both the MF Sim. and IRM suggests that the high frequency noise 
comes from the fine structure of large avalanches (T^ >> I//1, is the avalanche duration) 
and not small individual pulses (T^ ~ l/Zi)!^. Through simulation we verify this claim 
in the MF Sim. and IRM. When we subtract all avalanches smaller than T^ax/^: from the 
Barkhausen train and then determine the second spectra we find no change in V2. However, 
when we subtract all avalanches larger than from the Barkhausen train, we find that 

the second spectra flattens and that there is an evident separation between the 1.5 spectra 
curves (i.e. increased /i dependence). 

Our study of higher order spectra is a powerful tool to further our understanding of noise 
in disordered systems. Our mean field results are applicable to a large array of systems, in 
particular systems discussed in Q]. In a future long paper we plan to include details of all 
the calculations, and a comparison with experiment results. 
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FIG. 1: We present r2(A^) in the MF Sim. and in the IRM. For high frequency T2{N) 
Inset: Ti/2{N) for the MF and the IRM. For high frequency Ti/2{N) ~ see Tabled 
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FIG. 2: We present the -Re{S'i.5(/2, /i)} in the MF Sim. and in the IRM. At high frequency 
^e{S'i.5(/2, /i)} ~ f2^^'^ (go to Tabled). There is no visible flattening present due to r^^^(/i) 
{9 = 0) term, since the magnitude of this term is small relative to the /2 dependent term. 
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FIG. 3: We present the 82(12, fi) in the MF Sim. and IRM. Notice the flattening due to the 
^2^\fi) = 0) term. At high frequency 5'2(/2,/i) ~ /2^^^+ (DC background). We use a non- 
Unear curve fit of the form: AO * x""^^ + A2 to determine the V2 exponent given in Tabled Inset: 
5*2^ (/2, /i) is the second spectra in the MF Sim. and the IRM with the theoretical DC background 
term subtracted. The bold lines adjacent to the MF curve and the IRM curve have an exponent 
of -2 and -1.8, respectively. Please note that these curves (in the inset) do not show the full range 
of data points, and are simply presented to illustrate that we have a power law in the absence of 
the DC background term. 
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